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Abstract 

The standard analytic solution to the Dokshitzer-Gribov-Lipatov-Altarelh-Parisi (DGLAP) 
equation in Melhn space is improved by resumming the large x divergences. Explicit results are 
given to next-to-leading order and next-to-leading logarithmic accuracy. Numerically, the theoret- 
ical error is found to be reduced by the resummation for a large range of x. 
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Parton distribution functions (PDFs) are indispensable for the perturbative calculation of 
all high energy processes involving incoming hadrons. In particular, their precise knowledge 
is essential for the successful interpretation of the LHC experiments. Likewise, fragmentation 
functions (FFs) are required for inclusive hadron production calculations. By solving the 
spacelike/timelike DGLAP equation [l| 

j^D{N, Q') = PiN, a,iQ'))D{N, Q') (1) 

for the vector D of PDFs/FFs (the formulae of this letter apply to either causality, which 
is therefore not specified), their evolution is determined since the matrix of splitting 
functions P can be calculated perturbatively as P = J2'^=i^lP'^^~^^ y where = a<(/(27r), 
a. being the strong coupling. As P reaches next-next-to-leading order (NNLO) precision 
2j, resummation of divergences at phase space extremes becomes relevant. This letter is 
concerned with soft gluon divergences at large x, which have a general relevance since their 
effect propagates to all x through the evolution. For example, we show later that, relative 
to standard NLO evolution, the reduction in the theoretical uncertainty due to large x 
resummation is comparable to that on going to NNLO for a large range of x. In general, 
these divergences take the form a"[ln"~^(l — x)/{l ~ x)]^ (or a" In"'*'"^"'' in Mellin space), 
and r = 0, n labels the class of divergence. They factor out of a hard subprocess according 
to 3 

Wias,N) = W,U(^s,N) (^a:W^'SiN)j , (2) 

where the fixed order (FO) series, in parenthesis, is free of these divergences since they 
are all contained in W^es- At large A^, W is approximated by Wres if its divergences are 
resummed, which involves writing it as an exponential, then, in the exponent, grouping 
term, of the =an,e cla., giving a .e™= in a. keeping a.lniV fixed. Such expan.on. exi.t 
for coefficient functions for various processes, e.g. inclusive hadron production |J], deeply 
inelastic scattering and the Drell-Yan process j^. 

In this letter we present an analytic solution to the DGLAP equation which improves the 
FO accuracy by resumming large x divergences. We work in a minimal subtraction (MS) 
scheme where P is resummed: its diagonal components approach [1/(1 — x)]+ (InA^) while 
the off-diagonal components become constant [6]. Our approach is based on the commonly 
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used Mellin space solution The DGLAP evolution is written in the form 

D{N, Q2) = E{N, a,(Q2), as{Ql))D{N, Ql), (3) 

then the = ^^(Q^) and = as(Qg) dependences in the higher order parts of E are 
factored out, 

E{N, as, ao) = U{N, a,)^Lo(A^, "s, ao)U-\N, ao), (4) 



where i?LO is the LO evolution, given formally by 



Ei^o{N,as,ao) = exp 



(5) 



/3o ao 

Here, /3o is the first coefficient appearing in the perturbative expansion for the evolution of 
as, da,(/i2)/dln/i2 = -EZoPna^^+'i^')- Then 

— = -— + ^ rf/,p(°)l (6) 
das l3o f^otts 

replaces the DGLAP equation, where 

/ P p(o) \ 

(7) 



The accuracy of the evolution depends on how the matrix U is approximated. The standard 
approach is to use the the FO series 

oo 

UiN,as) = l + Y,U^^\N)a:, (8) 

n=l 

which can be determined order by order from Eq. ([6]) since the i?'-"^ depend on the 
for m < n — 1 only. However, the accuracy may be improved by resumming the large x 
divergences, which is the main task of this letter. The result is that at large A^, its diagonal 
components take the form Uii{N,as) — exp[^^j^ ^/j/ (a^ In A^)a^""'^] (recall that r, in the 
square brackets here, labels the class of divergence). 

The resummation of the non singlet (or valence quark) component of E is trivially 
implemented by solving the DGLAP equation exactly. However, we discuss it in detail 
since the resummation in the singlet evolution to be performed later is similar. Equa- 
tion (jlj) simplifies to the scalar equation i?Ns = {cLs/do) ^s/w (t4vfs(«s)/f^Ns(c^o))- In 
the unresummed case, U-^s{(^s)/U-^^s{(io) is expanded as a series in 0^,09, e.g. at NLO 
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it is taken as 1 — (a^ — aQ)R^l/ where R^l = Pj^g — (3q)P^^. In the re- 
summed case, it is calculated in the form given by Eq. ([2]), e.g. at NLO it becomes 
(^Ns('^s)/^Ns('^o)) 1 — (c^s — cio) (^f^NS^"* + ^Ns/ l-^oj ; whcrc U^s^'^ is the coefficient of the 
0{as) term in the FO expansion of f/j[||. To obtain the f/j[[g requires solving the equa- 
tion dU-^s/dcis = QnsU-ns at large N, where Qns = Pns/P + -Pns /(/^oOs). Of course this 
holds for all A^, as can be seen by substituting Eq. (jlj) into the DGLAP equation. For 
example, at NLO, the next-to-leading logarithms (NLLs), being r = 1 (0(as In A^)) terms, 
are resummed by taking f/j|| = exp[as[/Ng ^''] with ^/^g"^'' = —R^l/(3o + 0(1), and taking 
the last equation to be exact leads to the exact solution to the NLO DGLAP equation, 
^Ns = (a,/ao)-^NsV/3oexp[-(a, - ao)Ri^l/ (3o]. 

Analytic resummation of the singlet E is nontrivial, since its matrix structure rules out 
an exact and analytic solution to the DGLAP equation . We will present for the first time an 
analytic approach for resumming the FO singlet E. We first derive in detail the unresummed 
E, which will be modified later in order to implement the resummation. To put E'lo in Eq. 
(IS]) into a form which can be explicitly evaluated, we diagonalize P^^\ 

P(°) = A+M+ + A_M-, (9) 

where 
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p(0) , p(0) , J(p(.0) _ p(0)N2 _4p(0)p{0) 

gg ^ gg V ^9 ) ^^qg ^gq 

are the eigenvalues of P'^'^\ and the projection operators 



(10) 



M± = ^ [P(0) - A^l] (11) 



obey M±M^ = 0, M±M± = M±, ^.M^ = 1. Then 



Ei^o{N,as,ao) = Y,M\N)[^^ '° . (12) 

To diagonalize all occurrences of P'^^^ in [/, we work with the ijth. "components" projected 
out by operating on the left by A'P and on the right by MK Note that the sum of all 
"components" of any 2x2 matrix A gives back the full result, i.e. A = 'Y^^-M^AMK 
Matching coefficients of in Eq. ([6]) gives 

= V ^ — —M'R^"^M^. (13) 



\j - Ai - Pon 



(Note that P(o)M* = M'P(o) = \iM\ which follows from Eq. has been used here.) The 
right hand side is clearly a sum over the "components" M''U^"-^MK 
At NLO, 
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Aj — Aj — /3o 



where 



Pi 



^(1) ^ p(l) _ ^p(0)_ 

/3o 



(14) 



(15) 



In U , and therefore in general in there is a simple pole on the positive real axis of space 
when A+ — A_ — /9o = and 2 simple and 2 double poles off the real axis when A+ — A_ = 0, 
which will limit the choice of the contour used for the inverse Mellin transform. Fortunately, 
they manifest themselves in Eq. (jlj) as spurious NNLO terms, which may be subtracted. 
The NLO expansion in a<j, cto with Plo fixed. 



E = E^o + U^'^E^oas - E^oU^'^ 



ao, 



(16) 



is free of these poles, which is most easily seen in its projections: The "diagonal" {i = j) 
components read 



n \ ^0 
M'EM' = ( — 1 M' 



1 - (as - ao) 



PW 

(3o 



(17) 



and the "off-diagonal" {i 7^ j) components read 

1 



M'EM^ 



Aj — Aj — /3o 



(18) 



It is clear that the pole in U for which Aj — Aj — /5o = only appears in the "off-diagonal" 
component M^EM^ . Its cancellation can be seen by setting Aj = \i + Po + e and taking the 
limit e — > 0. The components of U containing simple and double poles for which A+ — A_ = 
appear in both the "off-diagonal" and "diagonal" components of E. The 1/(A+ — A_) terms 
cancel since E is invariant with respect to the interchange A_|_ ^ A_. The 1/(A_(_ — A_)^ 
terms, which give rise to both simple and double poles, can be seen to cancel by taking 
A+ = A- + e etc. as before. 

Now we incorporate the resummation in the solution. We require the asjTiiptotic be- 
haviour of U at large N. This is determined by using Eq. ([6]), or equivalently (fT3l) . to obtain 
the for all n as functions of the p(™'), m = 0, ...,n, whose asymptotic behaviour in 



MS schemes was given in the quark-gluon basis earher. Fortunately, this basis' components 
coincide with the projections we have been using: From Eq. ffTOj) . the large N behaviour of 
P implies that A+(„) ~ P'f] -.^ therefore 



This means that performing the projection MMM-' has the same effect as setting all com- 
ponents to zero except A^^ap where = <i{,g)- One immediately finds from Eq. f|T3|) 
that the "off-diagonal" components of U fall like 1/lnA^ and may therefore be left as a se- 
ries in ttsi while the "diagonal" components of U cannot be approximated in this way since 
M^U^"'^M^ grows like In" A^. To resum it, we solve the "diagonal" components of Eq. in 
the large N limit. The matrix multiplications on the right hand side can be decomposed 
according to AB = Y^--^{M^AM^){M^BM^), where B is also any 2x2 matrix. Then any 
product of an "off-diagonal" component of U with an "off-diagonal" component of R or P^°^ 
in either order contains no divergences and thus is neglected. Neglecting 0(1) terms, the 
result is 



where the matrix Qu = M' {P/(3 + P(°V(/5oa.)) at large N. This can be solved order by 
order to obtain a series for M^UM^ which can be identified and resummed. The diagonal 
component {Qn)aiai grows like InA^, while the other diagonal component falls like 1/lniV. 
The off-diagonal components both approach a constant. 

The NLO evolution contains divergences up to the NLL level. The r = (LL) terms are 
all contained in -Elo- The r = 1 terms are all contained in U, and according to Eq. fl2Ul) can 
be accounted for by taking 




(19) 



_d_ 
dtts 



{M'UM') = Qii{M'UM'), 



(20) 



M'UM' = Uf^M' + 0(a,(a, In AT)'"), 



(21) 



where the Ul are constrained at large by writing 



U 



exp agUl 



[i](i) 



(22) 



and requiring that 



[i](i) 



1 



(23) 
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The 0(1) terms in uj}^^^^ are not constrained, and will be chosen 



Ra}ai is calculated with Paia^ 



Pi%lnN in Eq. where 



ater. In this approximation, 



p{0) . 



(24) 



p{0,l) _'^^p(0,l) 
99 Ci. ■ 



From these results, and using the form in Eq. ([2]), the LL and NLL resummed U at NLO 
will be chosen as 



U = uf^M' 1 + a, (u^^^ - Uf-^^^h^ + a, ^ M'U'-^^M^. (25) 

The first line gives the resummed form for the "diagonal" components, while the "off- 
diagonal" components contained in the last line are as for no resummation. This agrees 
with the FO series when expanded to NLO, and gives the correct large behaviour, not- 
ing that the IniV divergence in M^U^^^'M^ cancels that in ul^^^^^M^. Analogously, is 
found to be equal to Eq. ( !25|) with f/^^-* and ul^^^^^ everywhere multiplied by —1. Then 
UU~^ = 1 + 0{a1) + 0(as(as In A^)™), as required. As in the unresummed case, the re- 
summed U and thus E contains problematic poles at A+ — A_ — /3o = and A+ — A_ = 0. 
To deal with this problem, and remain close to the unresummed approach of Ref. [tI], we 
expand E. Note however that uf^ must be kept fixed to preserve the resummation. Then 
the "diagonal" components read 



M'EM' 



1- {as- ao) 



(26) 



while the "off-diagonal" components read 



M'EM^ 



Aj — Aj — /3o 



ao 



ao 



'^Ul^ia. 



M'R^^^MK (27) 



These last two equations, together with Eq. (123!) . are our main results. They agree with 
Eq. f|T6|) when expanded to NLO without keeping uf^ fixed, i.e. when the resummation is 
undone, and obey the normalization condition 



E{N,as,as) = 1. 



(28) 



7 



The next step in eliminating these problematic poles involves tuning the f/]^"^^'' at finite 
A^. The infinite number of scheme choices refiects the uncertainty of the unknown higher 
order terms, a ubiquitous feature of perturbation theory. The pole cancellation is easier 
to achieve by imposing "common sense" constraints obeyed by the unresummed evolution, 
e.g. the requirements that ?7]^'*^^^ = at each pole so that the resummed evolution coincides 
with the unresummed one there, that the dependence of U occurs only through the X± 
and R^^\ and that the evolution is invariant under the interchange A+ ^ A_. Of course, 
there is also the resummation condition of Eq. fl23l) . One possibility that enforces all these 
conditions and leads to the pole cancellation is 

7-r[i](i) _ (A+ - A-)^ [(A+ - A-)^ - , . 

u^± -K\^ p^T+a:? ' ^ ^ 

where, for either i = ± (defining (C+,C_) = {Cf.Ca)), 

^=2^(^I^) (30) 
Finally, to study the effect of resummation on the theoretical error we vary the renor- 
malization scale /i by fixing /i^ = kQ^ and varying k. We also set /ip = kQl in oq as 
required by Eq. fl2Sl) . This means that 0^(0) = ^^(fcQ^o)) '^^^ expansion variables 

instead of a5(Q^g^), and (cis/ao)~'*''^^° and uf\as{())) are treated as fixed. The NLO relation 
<^s{Q\))) = as(o)(l + as(o)/?o In A;) implies the replacement P^^'' P^"'^^ + P'^°^/?o In A; everywhere. 
Note that P£l, -> - 2Cil3o In k in Eq. ([2 



For our numerical analysis, we fix the number of active quark fiavours nj = 5, put A^^^ = 
200 MeV and Qq = 2 GeV, and apply spacelike evolution to the non singlet distribution 
Dns{x,Qo) oc {1 — xY, being a typical large x behaviour. Resummation reduces the scale 
variation over a wide range of x (see Fig. [T]), as anticipated in the introduction. Away 
from the large x region, the scale variation is even reduced to around that of the NNLO 
calculation for low values in particular. For further analysis we refer the reader to a 
future publication where also resummation at NNLO will be studied. Fortunately, the three 
calculations accidentally give similar results around the commonly chosen value A; = 1, so 
that unresummed divergences should not have harmed global fits of PDFs/FFs. 

In conclusion, we have presented an analytic approach for resumming large x divergences 
in evolution. This makes a substantial numerical difference to the evolution away from 
k = 1. Since we perform the resummation in the "diagonalization" scheme of Ref. 
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FIG. 1: Renormalization scale variation of the non singlet evolved with our resummed NLO ap- 
proach (solid line), as well as with the usual unresummed NLO (dotted) and NNLO (dashed) 
approaches. 

that is used to solve the FO DGLAP equation, our approach is the simplest one. Since it 
provides both an important conceptual improvement to perturbative QCD and leads to a 
more accurate determination of PDFs and FFs, our approach is valuable for the description 
of phenomenology at the LHC and other coUiders of the foreseeable future. 

This work was supported in part by the German Federal Ministry for Education and 
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